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PREFACE 
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of  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under 
ARO  Projects  No.  PF403,  PF404,  and  P32S-22B.  The  authors 
of  this  report  were  K.  P.  Fewell  and  P.  A.  Kessel,  ARO,  Inc. 

The  data  analysis  was  completed  on  June  1,  1975,  and  the 
manuscript  (ARO  Control  No.  ARO-PWT-TR- 75-  107)  was  sub-  — 
mitted  for  publication  on  June  30,  1975. 
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1.0  INTRODUCTION 

In  tin-  past  several  years,  the  already  imposing  requ ireinonts  of 
tost  fueilities  designed  to  evaluate  the  ablation  response  of  materials 
have  been  farther  complicated  by  an  emerging  interest  in  erosion  by 
In  perveloeity  particle  impact.  Such  erosion  of  heat  shield  material  is 
known  to  occur  when  a reentry  vehicle  traverses  a cloud  of  solid  or 
liquid  particles  at  hypersonic  speeds.  These  clouds  may  be  of  natural 
or  man-made  origin  and  contain  particles  with  diameters  up  to  1,000  p , 
although  tiie  range  of  primary  interest  is  50  to  500  p . Hypervelocity 
impact  of  a stream  of  particles  on  a solid  surface  can  both  increase  the 
heat -transfer  rate  to  the  surface  and  remove  base  material  by  mechanical 
fracture  of  surface  elements. 


1 


The  basic  physical  processes  involved  in  particle  erosion  can  best 
be  studied  on  a single  impact  basis  in  a ballistic  range  facility.  In  such 
a facility,  the  velocity,  shape,  and  mass  of  the  impacting  body  can  be 
very  accurately  determined.  The  net  erosion  process,  however,  is 
composed  of  a great  number  of  successive  impacts  occurring  at  a high 
rate  on  a given  area  of  exposed  surface.  To  study  this  process,  test 
facilities  have  been  developed  which  produce  large  numbers  of  impacts. 
Particle  cloud  generators  have  been  developed  for  gun  ranges.  To  pro- 
duce erosion  simulation  for  times  longer  than  the  flight  time  in  a gun 
range,  particle  injection  systems  have  been  developed  for  use  in  wind 
tunnels.  In  the  AEDC/ PWT  Dust  Erosion  Tunnel  (DET),  for  example, 
solid  particles  of  50-  to  500 -p  diameter  are  accelerated  to  velocities 
of  4,000  to  (>,000  ft/sec.  To  attain  these  high  velocities  with  particles 
of  these  diameters,  it  is  necessary  to  drag-accelerate  them  in  a long, 
unconventional  nozzle  having  a very  slow  rate  of  expansion.  How- 
ever, even  though  it  is  an  arc-heated  facility,  the  DET  does  not 
qualify  as  an  ablation  facility  because  the  heating  rates  produced 
in  this  unconventional  nozzle  are  relatively  low. 


Since  erosion  and  ablation  often  occur  simultaneously  in  flight, 
there  is  considerable  interest  in  developing  a test  facility  having  a 
combined  capability.  Although  particle  velocities  of  18,000  ft/sec  are 
encountered  in  flight,  it  is  considered  that  velocities  of  8,000  ft/sec 
or  greater  would  be  of  interest  in  a ground  test  facility,  with  particle 
diameters  of  50-  to  500 -p  . As  an  outgrowth  of  the  RHEA  facility 
feasibility  studies  at  AEDC,  a facility  configuration  for  simultaneous 
ablation  and  erosion  testing  has  been  proposed,  and  this  report  describes 
the  analytical  and  experimental  investigations  for  evaluation  of  the  new 
concept. 
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Calculation  of  the  acceleration  of  particles  injected  directly  into 
the  nozzle  of  a conventional  arc-heated  ablation  facility  showed  that 
because  of  the  short  length  of  the  nozzle,  final  velocity  of  100-p  par- 
ticles would  be  less  than  one-half  the  gas  velocity  (Fig.  1).  The 
nozzles  for  typical  ablation  facilities  are  kept  to  a short  length  to  min- 
imize heat  losses  to  the  walls.  High  particle  velocity  and  high  heating 
rates  appear  to  be  mutually  incompatible  in  arc-heater  nozzles. 


Figure  1.  Gas  velocity,  particle  velocity,  and  temperature  in  an 
M = 2.0  ablation  nozzle,  D*  = 0.375  in. 


Studies  of  the  RHEA  facility  (Ref.  1),  an  advanced  concept  in- 
volving dynamic  interaction  of  particle  and  light-gas  flows,  focused 
attention  on  two  phenomena  which  have  since  been  incorporated  into  a ; 

proposed  ablation/erosion  facility.  First,  is  the  well-known  capacity  j 

of  a light  gas  to  provide  high  gas  velocities  for  relatively  moderate  | 

stagnation  temperatures.  Second,  as  indicated  by  general  particle 
trajectory  studies,  particles  on  the  order  of  100-p  diameter  traveling 
at  thousands  of  feet  per  second  are  highly  resistant  to  deflection  caused 
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by  aerodynamic  forces  imposed  by  most  Qow  fields  of  interest..  'That 
is,  tiie  trajectories  of  such  particles  w ill  be  very  nearly  linear.  A 
facility  concept  which  utilities  these  two  phenomena  is  shown  in  Fig.  2, 


In  this  proposed  facility,  ablation  capability  is  obtained  by  use  of 
a conventional  high-pressure  arc  heater  operating  with  air.  Particles 
are  accelerated  to  hypervelocities  in  a separate  light-gas  acceleration 
nozzle  and  then  transferred  to  the  ablation  jet  by  impinging  the  two 
flows  and  allowing  particle  momentum  to  effect  the  transfer  across  the 
mixing  layer.  Analysis  and  several  experimental  efforts  related  to  the 
feasibility  of  this  concept  are  reported  herein. 


2.0  COMPUTER  PROGRAMS  FOR  ANALYSIS  OF  GAS/PARTICLE  FLOWS 


Several  computer  programs  were  used  in  the  analysis  of  a combined 
ablation/ erosion  nozzle.  In  each  of  these,  the  particle  drag  coefficient 
is  an  important  parameter,  and  two  different  drag  correlations  were  used. 
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2.1  ONE-DIMENSIONAL  GAS/PARTICLE  FLOW  WITH 
THICK  BOUNDARY  LAYERS 


The  simplest  conceptual  approach  to  producing  high  speed  in  par- 
ticles entrained  in  a gas  is  to  inject  the  particles  upstream  of  the  throat 
of  a supersonic  nozzle  and  allow  drag  c.cce] eration  to  the  desired  velocity 
level.  A quasi-one -dimensional  computer  program  was  previously  devel- 
oped at  AEDC/VKF  for  calculating  such  particle  acceleration  in  supersonic 
nozzles  having  turbulent  boundary  layers  on  the  nozzle  wall  (Ref.  2).  This 
program  was  used  to  analyze  the  nozzle-particle  flow  shown  in  Fig.  1 for 
a conventional  low  Mach  number  nozzle  and  also  to  design  the  acceleration 
nozzle  of  the  DET.  In  the  latter  case,  it  was  found  that  optimization  of  a 
nozzle  for  drag  acceleration  of  particles  in  the  100-P  range  of  diameters 
requires  the  use  of  a very  long  nozzle,  which  results  in  very  thick  (nearly 
merged)  boundary  layers.  In  the  program,  the  boundary  layer  calcula- 
tion is  coupled  to  the  one -dimensional  channel  flow  calculation,  but  the 
particle  acceleration  calculation  is  not  coupled;  i.  e.  , the  particle 
acceleration  occurs  in  a predetermined  gas  dynamic  flow  field  and  does 
not  affect  that  flow  field. 


A solution  for  one -dimensional  channel  flow  with  wall  friction  and 
heat  transfer  is  obtained  by  integration  of  the  following  three  conservation 
equations  and  an  equation  of  state; 

Conservation  of  mass: 


*2  + dV  + dA  „ o 
p V A 


Conservation  of  momentum: 


;vfita=-2cfpt 


dx  dx 


D 


Conservation  of  energy: 


V ^ = . 4 2ll  (h  - h ) 

JJ j..  n " 


dx  dx 


Equation  of  state: 


H » h(p,p) 


(1) 


(2) 


(3) 


(4) 
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where  the  skin  friction  coefficient  and  the  Stanton  number  are  given  by  j 

the  foil  wing  relations  for  turbulent  flow:  j 


Cf 


Ch  = 0.53  Cf  (6) 

The  original  VKF  program  was  formulated  for  use  with  either  air, 
helium,  or  hydrogen  as  the  gas  flow.  The  properties  of  high  temper- 
ature air  are  obtained  from  a real  gas  subroutine  which  computes 
equilibrium  composition  and  thermodynamic  properties  with  correc- 
tions for  inter  molecular  forces.  Helium  is  treated  as  an  ideal  gas, 
and  hydrogen  is  modeled  as  a gas  which  vibrates  but  does  not  disso- 
ciate. The  constants  in  the  skin  friction  relation  were  determined 
from  an  empirical  correlation  of  test  data  on  turbulent,  flat  plate 
boundary  layers  and  are  given  in  Table  1. 


Table  1.  Constants  in  Skin  Friction  Relation  Eq.  (5) 


Integration  of  the  foregoing,  equations  gives  the  flow  solution  in 
terms  of  average  values  (p,  h,  V)  in  a channel  of  known  geometry, 
with  skin  friction  and  heat  transfer  at  the  walls  of  the  channel. 
Corresponding  "core  flow"  values  are  obtained  from  an  isentropic 
solution  (no  skin  friction  or  heat  transfer)  for  the  same  pressure 
distribution.  Then,  the  boundary -layer  growth  from  the  wall  to  the 
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core  is  obtained  by  integration  of  a differentia],  equation  for  momen- 
tum thickness  and  by  reference  to  two  supplemental  equations  for 
total  thickness  and  displacement  thickness: 

d0  f 1 dp  ^ 6*,  1 dV  1 dDl  Cf 

— + 0 + (2  + — ) + = — 

dx  p dx  9 V dx  D dx  2 


(1  _ 1$1)2  = £_v 


6 = 6*  + 6.4  8 (c 

The  relation  for  total  thickness  was  based  on  experimental  data  in 
long,  narrow  channels.  The  complete  quasi-one-dimensional  solution 
for  the  channel  is  obtained  by  solution  of  Eqs.  (1)  to  (9).  The  solution 
is  valid  until  the  boundary -layer  growth  is  great  enough  that  the  total 
thickness  merges  on  the  flow  centerline.  The  increase  in  entropy 
which  would  occur  downstream  of  this  point  is  not  allowed  for  in  the 
program  as  presently  constructed. 

The  behavior  of  particles  in  the  given  flow  field  is  obtained  by 
integration  of  differential  equations  for  acceleration  and  particle 
temperature: 

dVp  _ Cp  AP  p(V-V„>  X [ <V-y|  (10 

dx  2 Vp  nip 


dTn  Q 

P c 


dx  (pv  C)  V 
P p P 


P(V-Vp)  StT°-V  ChP  Sp  (11, 

(pv  Cp)p  vp 


where 


Q a ^NUj  — ^ ^ * Q • 6 
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by  assuming  Pr 


P(V- v )0p 

Re  = — u: 

* y 

and  the  gas  temperature  is  based  on  the  slip  velocity,  The  particle 
Nusselt  number  was  obtained  from  Ref.  3 and  is  strictly  valid  only 
for  subsonic  slip  Mach  numbers.  The  relatively  large  p:  rticle  diam- 
eters used  in  the  present  study  result  in  supersonic  slip  velocities. 

To  allow  for  this,  the  particle  Reynolds  number  and  heat  transfer 
are  simply  evaluated  using  total  conditions  behind  the  particle  bow 
shock. 


The  program  version  used  in  the  present  study  was  modified  from 
that  of  Ref.  2 by  the  addition  of  the  particle  temperature  history  (Eq. 
(11))  and  by  extension  of  the  particle  acceleration  and  temperature 
solution  to  an  arbitrary  point  in  the  subsonic  region  upstream  of  the 
nozzle  throat.  The  boundary -layer  growth,  however,  is  still  assumed 
to  begin  at  the  nozzle  throat  m order  to  simplify  prediction  of  mass 
flow  in  the  nozzle. 


SKs***4* 


2.2  TWO-DIMENSIONAL  TRAJECTORY  CALCULATIONS 
IN  PARTICLE/GAS  FLOW 

In  the  analysis  of  feasibility  of  a combined  ablation  and  erosion 
facility,  it  was  necessary  to  perform  particle/flow  calculations  for  the 
case  in  which  the  slip  velocity  is  not  aligned  with  the  particle  velocity. 
In  this  case,  the  particles  follow  cmved  paths,  arid  a two-dimensional 
program  is  required.  In  this  program,  as  in  the  one -dimensional  pro- 
gram, the  particle  accelerations  are  computed  in  a predetermined  flow 
field,  with  no  coupling  of  the  particle  motion  upon  the  flow  field.  The 
particle  trajectories  are  obtained  by  integration  of  the  vector  equation 
of  particle  acceleration: 


= ,2  SE  JE-  v.»)  i 

dt  4 Dp  Pp  vs 

where  Vs  is  the  slip  velocity  vector  - Vp).  Since  the  drag  force 
and  acceleration  vector  are  in  the  direction  of  Vs,  the  components  of 
acceleration  are: 

dVPx  Vsx  dVp 

— ^ = (_JL)  _ £ (1 
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In  this  integration,  it  is  considered  that  Dp  and  p are  constant,  but  that 
Cj),  p,  and  Vg  are  functions  of  position  coordinates  depending  on  the 
nature  of  the  flow  field.  Programs  for  two-dimensional  trajectory  cal- 
culations have  been  constructed  for  the  following  flow  fields: 

a.  Uniform  flow. 

b.  Prandtl- Meyer  expansion  flow. 

Conditions  along  successive  expansion  rays  are  calculated 
by  mtegrating  the  following  expression  from  known  initial 
conditions: 


= exp  cos 


(— ) - cor 


-■<i>  / =&= 

M,  '/mTTT 


where  a is  the  angle  of  an  expansion  ray  measured  from  the 
horizontal.  The  integration  is  performed  numerically,  since 
the  real  gas  effects  which  dominate  the  area  of  interest  pre- 
clude analytical  integration. 

c.  Ablation/ erosion  nozzle  flow. 

In  one  comprehensive  program,  a real  gas  isentropic  expansion 
of  air  from  high  temperatures  is  calculated  and  matched  in  static 
pressure  to  a helium  or  hydrogen  flov'  downstream  of  an  oblique 
shock.  Then,  in  the  same  program,  conditions  upstream  of  the 
shock  in  the  light  gas  are  calculated  and  also  conditions  down- 
stream of  a similar  shock  in  the  airflow  are  calculated.  Finally, 
trajectories  of  particles  having  different  velocities  and  diameters 
uie  calculated  for  traversal  of  the  four  separate  but  adjoining 
flow  fie.1  as. 

2,3  PARTICLE  DRAG  COEir FtC!ENTS 

For  lack  of  a more  workable  procedure,  the  drag  coefficients  used 
in  most  analyses  of  particle  flows  are  assumed  to  be  those  of  smooth 
spheres.  In  many  cases,  the  particle  diameter  is  such  that  the  slip  veloc- 
ities are  low  subsonic,  in  which  case  the  well-established  incompressible 
sphere  drag  curve  giving  Cjj  as  r:  function  only  of  Re^  is  used.  This  curve 
becomes  asymptotic  to  the  Stokes  drag  coefficient,  = 24/Rej-j,  at  very 
low  Reynolds  numbers  «'Fig.  3).  In  t.he  present  case,  the  particle  diame- 
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Figure  3.  Sphere  drag  coefficient  versus  Reynolds  number. 


ters  were  such  that  supersonic  slip  velocities  were  often  encountered, 
and  it  was  necessary  to  recognize  the  Mach  number  dependence  as  well 
as  the  Reynolds  number  dependence  of  the  drag  coefficient. 

2.3.1  Crowe  Drag  Correlation 

One  widely  used  relation  for  sphere  drag  coefficient  in  the  compress- 
ible case  is  that  published  by  Crowe  (Ref.  4).  This  correlation  was  based 
on  sphere  drag  data  from  large  number  of  different  test  facilities  and  is 
as  follows: 


CD  " (cDinc  - 2)  exP  "3.07  f 0 • 5 ~ g (Re) 


logi  o g (Re) 


= 1,25  [ 


1 + tanh(0.77  logi  0 Re 


- 1.92)] 


(18) 


(17) 

r.  §®i 

+ 2 

L 2MJ 

1 
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where  CDinc  is  the  (nearly)  classical  incompressible  drag  coefficient  of 
a sphere.  Most  of  the  preliminary  analyses  of  this  report  utilized  this 
drag  correlation.  It  was  found,  however,  that  in  some  regimes  the  Crowe 
relation  did  not  agree  with  AEDC  sphere  drag  measurements  (Fig.  3). 

2.3.2  Correlation  of  AEDC  Sphere  Drag  Data 

A large  amount  of  sphere  drag  test  data  has  been  obtained  at  AEDC 
during  the  last  decade,  covering  a very  large  range  in  Mach  number  and 
Reynolds  number  in  the  same  test  facility,  a fact  which  results  in  a high 
degree  of  internal  consistency  and,  inferentially,  a high  level  of  absolute 
accuracy.  In  Ref.  5,  data  were  published  for  0.  1 < M < 6.  0 and  20  < Re^ 

< 100,000,  with  effects  of  sphere  temperature  given  for  M = 2 to  6. 
Further  data  were  reported  in  Ref.  6 for  0.  052  < M < 0.  80  and  0.  03  < 

Rejj  < 1000.  A comparison  of  these  data  (supplemented  by  Ref.  7)  with 
Eq.  (17)  evaluated  for  M = 1 is  shown  in  Fig.  3.  As  indicated  here  and 
in  Ref.  9 as  well,  the  use  of  Eq.  (17)  can  result  in  drag  coefficients  which 
differ  substantially  from  test  data  at  Rej-j  below  100. 

A curve  fit  was,  therefore,  obtained  for  the  AEDC  test  data: 

CD  = exp  | a2(8n  ReD)  2 + ai  Cn  ReD  + aQ|  (1  + k 8n  -£■)  (20) 

with  different  coefficients  obtained  for  various  Mach  number  ranges  (Table 
2)  and  where  k is  a function  of  both  M and  Re^. 

Table  2.  Interim  Curve-Fit  to  Sphere  Drag  Coefficients 


(VKF  RANGE  DATA,  REF.  6)  3?  • bjM?  4 bjM  4 b0 

M ■ 0. 75-4.2  Re0  ■ 20  !o  5000  cQ  ■ exp  ja;>Un  Rec>2  ♦ 3[/n  ReD  + a0j  aj  • 03 M2  ♦ CjM  4 c0 

a0  ■ d?M2  + djM  * dg 


M 
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A comparison  of  the  curve  fit  with  test  data  is  given  in  Fig.  4 for  two 
different  Mach  numbers. 


Figure  4.  Comparison  of  curve  fit  and  AEDC  data  on  sphere  drag  coefficient. 


2.3.3  Nature  of  "Real"  Particles 

Both  of  the  empirical  drag  correlations  described  above  are  strictly 
valid  only  for  single,  smooth,  spherical  particles.  However,  real  parti- 
cles used  in  testing  often  display  very  nonspherical  shapes  and  a fairly 
broad  distribution  of  sizes  (Fig.  5 and  G).  In  practice,  then,  one  would 
expect  to  observe  a dispersion  of  results  about  a "mean"  velocity  or 
trajectory  calculated  for  an  effective  average  spherical  particle.  Note 
that  it  would  be  possible  for  effects  of  nonsphericity  on  drag  to  be  com- 
pensated by  errors  in  determination  of  the  particle  sizes. 

The  particle  samples  formed  of  small  glass  beads  (Figs.  5 and  6) 
were  used  to  reduce  the  dispersion  of  results  for  careful  calibration 
work.  Experimental  measurements  of  particle  behavior  invariably  show 
this  reduction  in  dispersion. 


3.0  EXPERIMENTAL  EVALUATION  OF  GAS/PARTICLE 
FLOW  PROGRAMS 


Adequacy  of  the  gas/ particle  flow  programs  as  design  tools  was  H 

evaluated  in  several  different  experimental  programs,  some  designed  j 

specifically  for  the  evaluation  and  some  designed  for  other  purposes.  ! 
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Figure  5.  Photomicrographs  of  particle  samples  showing  individual  shapes. 


3.1  EVALUATION  OF  THE  ONE  DIMENSIONAL  PARTICLE 
ACCELERATION  PROGRAM 

Uncertainties  contained  in  the  particle  acceleration  program  in- 
clude the  drag  coefficient  correlation,  the  boundary -layer  calculation 
for  very  thick  boundary  layers,  and  in  the  case  of  air,  the  high  tem- 
perature thermodynamic  property  routine.  The  experimental  evaluations 
were  made  under  conditions  in  which  the  various  uncertainties  were 
present  at  the  same  time.  Therefore,  only  ar^  overall  evaluation  of  the 
program  was  possible,  not  a separate  evaluation  of  the  individual 
uncertainties. 
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a.  100-g  MqO  particles  and  100-pt  glass  beads  j 

Figure  G.  Particle  size  histograms.  4 


3.1.1  PWT  Dust  Erosion  Tunnel  (DET) 


The  DET  was  designed  in  1970  as  an  arc-heated  wind  tunnel  having  i 
high-speed  particle  capability.  The  acceleration  nozzle  for  this  * 

facility  was  designed  using  the  unmodified  version  of  the  one- dimensional  | 
program  (Section  2.  l).  This  nozzle  has  a throat  diameter  of  0.  5ti3  in.  , | 

ah  exit  diameier  of  15.  34  in.  , and  length  of  209  in.  The  modified  pro-  1 

gram  was  recently  applied  to  this  nozzle,  and  the  results  are  given  in  J 

Fig.  7. 


AE DC  TR *75-1  54 


[ 

! 

[ 


I 

i 

i 

; 


i 

j 

i 

j 

fj 

£ 


b.  1-fi  Ti02  particles 
Figure  6.  Concluded. 


The  gas  velocity,  the  velocity  and  temperature  of  100-p  magnesium  oxide 
particles,  and  the  gas  Mach  number  are  shown  as  functions  of  distance 
along  the  nozzle  for  pQ  = 1,000  psi  and  HQ  = 1,446  Btu/lbm, 

The  experimental  Mach  number  at  the  nozzle  exit  was  determined 
by  measurement  of  the  impact  pressure  and  the  nozzle  static  pressure 
and  was  found  to  be  very  close  to  the  value  predicted  by  the  one-dimen- 
sional program  (Fig.  7b).  The  impact  pressure  profiles  (not  shown 
her*'  also  verify  that  the  exit  flow  is  very  close  to  the  merged  boundary- 
layer  condition. 

Particle  velocities  were  measured  at  the  nozzle  exit  by  double - 
pulse  optical  holography,  and  the  dispersion  in  velocity  because  of  varia- 
tions in  particle  size  was  clearly  observed.  However,  the  mean  value 
of  4,  900  ft/ sec  was  very  close  to  the  value  of  5,  000  ft/ sec  predicted 
for  100-p  particles.  The  predicted  gas  velocity  was  8,  250  ft/sec.  The 
velocities  measured  for  different  diameter  particles  are  compared 
with  the  one-dimensional  program  predictions  using  both  drag  coef- 
ficient correlations  in  Fig.  8. 
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CONTOUR  OF  SUPERSONIC  SECTION  OF  DET  NOZZLE- 
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b.  Supersonic  section 

Figure  7.  Gas  velocity,  Mach  number,  particle  velocity,  and 
temperature  in  DET  nozzle,  D*  = 0.563  in. 
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Figure  8.  Comparison  of  calculated  and  measured  particle  velocities  in  DET. 


3.1.2  Light-Gas  Acceleration  Nozzle 

To  evaluate  the  performance  of  the  light -gas  accelerator  in  the 
proposed  ablation/ erosion  nozzle,  the  one -dimensional  program  was 
applied  to  a small  nozzle  having  a large  L/D  ratio.  This  nozzle  was  a 
simple  conical  nozzle  having  a 0.  5-deg  divergence  half-angle  and  a 
throat  diameter  of  0.  119  in.  Results  of  the  calculation  for  a helium 
expansion  from  pQ  = 1,000  psi  and  TQ  = 2,000°R  are  shown  in  Fig.  9, 
for  a range  of  particle  diameters  from  10  to  1,000  M . The  nozzle 
dimensions  and  reservoir  conditions  for  this  calculation  were  selected 
on  the  basis  of  equipment  availability.  For  100-p  particles,  a velocity 
of  4,  500  ft/sec  is  attained  at  a length  of  18  in. , very  close  to  that 
attained  in  the  DET  in  a length  of  over  200  in.  and  for  an  accelerating 
gas  stagnation  temperature  only  40  percent  that  of  the  DET.  The 
calculation  also  indicates,  however,  a boundary -layer  merge  at  a 
distance  of  14  in.  As  noted  in  Section  2.  1,  the  gas  dynamic  part  of 
the  program  is  not  strictly  valid  downstream  of  this  point;  however, 
it  is  believed  that  the  particle  acceleration  calculations  will  not  be 
seriously  affected  for  an  appreciable  distance  downstream  of  this  point. 

In  Fig.  10  is  shown  the  particle  velocity  variation  in  the  same 
nozzle,  calculated  with  both  of  the  drag  coefficient  correlations  of 
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Section  2.3.  The  drag  correlation  based  on  the  AEDC  test  data  (Eq. 

(20))  gives  particle  velocities  approximately  10  percent  higher  than 
does  the  correlation  of  Eq.  (17). 

3.1. 2.1  Experimental  Apparatus  for  Light-Gas  Acceleration 

An  experimental  apparatus  was  assembled,  including  an  available 
electric  resistance  heater,  a stilling  chamber  and  particle  injection 
system,  a particle  acceleration  nozzle,  various  instrumentation  systems, 
and  a helium  and  air  gas  supply  system. 

The  resistance  heater  was  essentially  a carbon  steel  shell  with  sin 
inner  stainless  steel  liner  in  which  a spiral -wound  coil  of  Kanthal®  wire 
was  mounted,  properly  insulated  with  Fibre  Frax®.  The  design  condi- 
tions for  this  heater  were  1,000  psi  and  2,  000°R.  Electrical  resistance 
of  the  heating  element  was  approximately  2 ohms.  The  heater  was 
operated  with  a d-c  power  supply  having  two  ranges  of  operation:  4,000 
amp  at  250  v and  2,  000  amp  at  500  v.  A ballast  resistor  of  0.  4 ohm  was 
used  in  the  heater  circuit.  Performance  of  this  heater  was  calibrated 
with  air  for  mass  flows  of  0.  16  to  0.  3 lbm/sec  and  output  temperatures 
up  to  1,  800°R.  The  heater  was  not  calibrated  with  helium  because  of  the 
cost  involved,  but  results  observed  in  operation  of  the  acceleration  nozzle 
with  helium  indicated  that  approximately  the  same  temperatures  were 
attained  as  with  air. 

The  l-in.-ID  stilling  chamber  (Fig.  11)  was  machined  from  a stain- 
less steel  bar  and  was  attached  to  the  resistance  heater  by  a short  length 
of  stainless  steel  tubing.  Provision  is  made  in  the  stilling  chamber  for 


Figure  11.  Uoss  section  of  stilling  chamber  and  nozzle  of 
light-gas  accelerator. 
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measurement  of  stagnation  pressure  and  temperature  and  for  introduction 
of  particles. 

The  particle  acceleration  nozzle  itself  (Fig.  11)  was  made  from  a 
tapered  tube  welded  to  a throat  assembly.  The  tapered  tubes  obtained 
were  of  Type  347  stainless  steel,  having  a wall  thickness  of  0.035  in. 
and  a length  of  approximately  25  in.  The  tube  was  cut  to  the  required 
size  and  welded  to  the  pre-machined  throat  block.  The  nozzle  inlet  is  a 
circular  arc  transition  to  the  0.  5 -deg  divergence  angle  at  the  tapered 
tube.  The  throat  diameter  is  0.  120  in.  and  the  diameter  at  the  25-in. 
length  is  0.55  in. 

A simple  particle  mjection  system  was  constructed,  having  a 
carbon- steel  reservoir  with  an  inverted -cone  lower  surface.  This  tank 
was  separated  from  the  stilling  chamber  by  a ball  valve  and  metering 
orifice,  pressurized  by  a regulated  helium  supply.  Particle  injection 
is  90  deg  to  the  axis  of  the  flow;  no  attempt  was  made  to  inject  with  a 
velocity  component  parrallel  to  the  flow.  A vibrator  was  added  to  the 
system  to  ensure  a smooth  and  continuous  particle  flow.  The  initial 
configuration  of  this  system  provided  excessive  particle  flow  rates, 
but  after  modification  to  allow  injection  under  essentially  zero  pressure 
difference  acrors  the  orifices,  flow  rates  as  low  as  0.0006  lbm/sec 
were  obtained. 


Provisions  were  made  for  measuring  the  stagnation  conditions 
in  the  reservoir.  To  calibrate  the  exit  flow,  a three-probe  rake  was 
designedfor  traversal  of  the  flow.  This  rake  included  a metal -sheathed, 
Chromer -Alumel®  stagnation  point  temperature  probe,  a conventional 
impact  pressure  probe,  and  a static  pressure  probe.  The  static  pressure 
probe  was  eventually  removed  since  very  little  data  were  ever  retrieved 
from  it.  The  rake  was  driven  by  a low -speed  traversing  system. 

Photographs  of  various  components  of  the  installation  are  given  in 
Fig.  12. 

3. 1.2. 2 Aerodynamic  Calibration 

The  light-gas  acceleration  nozzie  was  calibrated  with  both  air  and 
helium.  Changes  in  nozzle  length  due  to  thermal  expansion  required  that 
the  calibration  probes  be  pre -positioned  some  distance  downstream  of 
the  nozzle  exit.  After  a 4 -min  stabilization  period,  the  gap  between  the 
probes  and  nozzle  was  narrowed  to  the  desired  value  of  0.  10  in. 
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a.  Heater  end  of  apparatus 


b.  Exit  of  nozzle 

Figure  12.  Photographs  of  light-gas  particle  accelerator  apparatus. 
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A typical  impact  pressure  profile  measured  at  the  nozzle  exit  for 
air  at  pQ  = 935  psi  and  T0  = i,  284°R  is  given  in  Fig.  13.  At  this  lower 


Tq,  the  nozzle  flow  should  not  be  merged  and  a small  region  of  constant 


impact  pressure  does  appear  at  the  centerline  of  the  profile.  Under  this 
condition,  it  is  legitimate  to  infer  a nozzle  exit  static  pressure  from  the 
impact  pressure  and  an  assumed  isentropic  expansion  on  the  nozzle 
centerline.  This  pressure  was  assumed  to  be  (essentially)  constant  over 


the  profile  ana  p/  p^  used  to  obtain  the  Mach  number  profile,  also  shown 


in  Fig.  13.  The  oscillation  in  the  impact  pressure  profile  is  believed 
to  be  the  result  of  oblique  shocks  at  the  nozzle  exit  due  to  overexpanded 
flow  conditions. 
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Figure  13.  Pitot  pressure  and  Mach  number  profiles  in  light-gas  particle 
acceleration  nozzle  operated  with  air. 


At  all  lower  stagnation  temperatures,  the  boundary  layer  will  be 
even  thinner,  so  that  the  same  procedure  should  be  valid  to  obtain  the 


centerline  Mach  number,  shown  in  Fig.  14  for  air  as  a function  of  TQ. 


The  experimental  M^  curve  is  parallel  to  the  inviscid  curve  (which 
decreases  with  TQ  because  of  real  gas  effects  in  air)  and  roughly  half- 
way between  the  inviscid  curve  and  the  viscous  flow  solution  obtained 
with  the  one -dimensional  program.  The  air  liquefaction  boundary  is 
shown  for  reference  as  well  as  velocities  calculated  for  50-  and  100 -d 
particles  of  MgO. 
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Figure  14.  Mach  number  and  calculated  particle  velocity  in  light-gas 
accelerator,  air. 

Corresponding  aerodynamic  calibration  data  obtained  with  helium 
at  pQ  = 750  psi  are  given  in  Fig.  15.  The  experimental  Mach  numbers 
again  fall  between  the  inviscid  line  (constant  for  a perfect  gas)  and  the 
viscous  nozzle  calculation;  however,  in  the  case  of  helium,  the  results 
are  closer  to  the  viscous  calculation.  The  trends  are  the  same  as  for 
the  calibration  with  air. 

3. 1.2.3  Particle  Velocity  Calibration 

During  the  initial  work  with  the  light -gas  acceleration  nozzle,  the 
only  technique  available  to  obtain  an  indication  of  high  velocities  was  that 
of  impact  specimens.  Photographs  of  two  1/ 16-in.  -diam  steel  rods  ex- 
posed to  particle  flows  for  periods  of  1 sec  are  shown  in  Fig.  16.  Note 
the  difference  between  the  specimen  exposed  to  2,  700-ft/sec  particles  of 
MgO  accelerated  by  air  and  the  specimen  exposed  to  5,  200-ft/sec  parti- 
cles of  MgO  accelerated  by  helium.  Even  specimens  exposed  to  100-p 
glass  beads  accelerated  by  cold  helium  display  characteristic  hyperve- 
locity impact  craters  (Fig.  17). 

An  extensive  effort  was  made  to  obtain  particle  velocity  measure- 
ments by  use  of  the  double-pulse  holography  technique  used  so  success- 
fully in  the  DET.  However,  the  much  higher  density  and  turbulence 


28 


AEDC.tr  75-154 


Figure  15.  Mach  number  and  calculated  particle  velocity  in  light-gas 
accelerator,  helium. 


Figure  16.  Photographs  of  impact  bars  after  1-sec  exposure  to 
50-p  MgO  particles,  air  and  helium. 
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Figure  17.  Photograph  of  impact  bar  after  exposure  to  100-m  glass  beads 
at  3,200  ft/sec  (calculated). 

levels  in  the  light-gas  particle  acceleration  nozzle  were  found  to  degrade 
the  object  beam  to  such  an  extent  that  the  holograms  were  not  recon- 
structive. Near  the  end  of  the  project  effort,  however,  velocity  meas- 
urements were  obtained  with  a laser  (doppler)  velocimeter  which  was 
optimized  for  the  conditions  of  the  light-gas  nozzle.  As  noted  below, 
the  velocities  measured  in  the  nozzle  were  close  to  predicted  values. 

3. 1.2.4  Evaluation  of  One-Dimensional  Program 

Ultimately,  the  one -dimensional  program  should  be  judged  on  the 
basis  of  its  intended  use,  namely,  to  predict  the  velocities  of  particles 
which  are  accelerated  by  drag  forces  in  special  acceleration  nozzles. 

A summary  plot  of  particle  velocities  measured  by  two  different  tech- 
niques in  two  different  facilities,  in  comparison  with  values  predicted 
by  the  one -dimensioned,  program,  is  given  in  Fig.  18.  This  comparison 
clearly  establishes  the  one -dimensional  program  as  a valuable  design 
tool  and  appears  to  eliminate  the  possibility  of  any  major  errors  in  the 
drag  coefficient  routine,  the  thick  boundary -layer  formulation,  or  the 
real  gas  routines  for  air  properties.  In  addition,  the  agreement  of  the 
aerodynamic  calibration  with  predictions  of  the  program  is  further  evi- 
dence that  the  gas  dynamic  formulation  of  the  program  is  sound. 
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Figure  18.  Summary  of  particle  velocity  measurements  in  DET  and  light-gas 
particle  acceleration  nozzle,  100-m  MgO  particles. 


3.2  7WO-DIMENSIONSL  TRAJECTORY  EXPERIMENT  IN  THE 
DUST  EROSION  TUNNEL 

Analytical  trajectories  computed  for  particles  in  various  flow 
fields  suggested  that  particles  of  100 diameter,  when  moving  at 
velocities  of  thousands  of  feet  per  second,  are  very  insensitive  to 
deflection  by  the  flow  fields.  Since  such  behavior  could  be  of  very 
great  use  in  operation  of  a combined  ablation/ erosion  nozzle,  exper- 
imental verification  of  this  was  considered  to  be  of  substantial  importance. 
An  experiment  was  designed  for  this  purpose  in  which  particles  would  be 
subjected  to  a cross  flow  created  in  Prandtl- Meyer  expansion  formed  in 
the  DET. 
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The  DET  facility  (Fig,  19)  provides  particle  flow  at  high  velocity 
in  a nozzle  which  may  be  operated  at  several  lengths,  the  extremes 
producting  Mach  numbers  of  5.  9 and  9.47.  This  facility  operates  with 
a 5 -MW  arc  heater  providing  reservoir  conditions  of  1,000  psi  and 
enthalpy  from  400  to  2,000  Btu/lbm. 
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Figure  19.  Schematic  of  AEDC/PWT  Dust  Erosion  Tunnel  (DET). 

An  adjustable  wedge  model  was  designed  and  fabricated  from 
stainless  steel  (Fig.  20).  By  varying  the  angle  of  attack  of  this  model, 
a Pr and tl- Meyer  expansion  flow  of  desired  strength  could  be  generated 
on  the  lee  side  of  the  model.  The  wedge  tip  is  replaceable,  and  up  to 
ten  3/ 16 -in.  -diam  impact  bars  were  mounted  on  the  lee  surface.  The 
trajectory  of  particles  in  the  Prandtl- Meyer  flow  region  can  be  deter- 
mined by  location  of  the  impact  boundaries  on  each  successive  bar.  To 
the  extent  that  the  particles  w 11  have  a certain  spread  in  size  (Fig.  6), 
shape  (Fig.  5),  and  probably  small  differences  in  initial  velocity  vector, 
the  boundary  of  the  impact  zone  will  be  somewhat  diffused,  increasingly 
so  at  the  more  downstream  impact  bars. 
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Figure  20.  Photograph  of  Prandtl-Meyer  expansion  test  model. 
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Twelve  tests  of  the  Prandtl-  Meyer  flow  model  were  run  in  the  DET 
(Table  3).  The  tests  performed  at  M - 7.  5 were  parasite  tests  performed 
during  operation  of  the  DET  for  another  purpose.  They  were  performed 
for  familiarization  only  since  prior  calculations  indicated  the  particle 
deflections  at  the  M-  7.  5 condition  would  not  be  detectable.  When  the 
DET  became  available  for  a test  dedicated  to  the  particle  trajectory 
study,  the  model  was  moved  upstream  to  the  M = 5.  5 station,  for  which 
conditions  the  theoretical  analysis  indicated  that  detectable  deflections 
over  a 7-in.  length  would  occur.  As  shown  in  Table  3,  six  tests  were 
performed  at  M = 5.  5,  using  glass  beads  of  two  different  diameters, 

MgO  particles,  and  Ti02  particles. 


Table  3.  Conditions  for  Prandti-Meyer  Flow  Tests  in  Dust  Erosion  Tunnel 


A photograph  of  the  model  after  a run  is  shown  in  Fig.  21  and  a 
photomicrograph  of  one  of  the  forward  impact  bars  is  given  in  Fig.  22. 
The  edge  of  the  impact  area  is  clearly  evident,  even  though  it  is  not 
perfectly  sharp.  The  trajectory  was  defined  by  the  mean  impact  point, 
assumed  to  lie  halfway  between  the  point  where  there  are  essentially  no 
impacts  and  the  point  where  there  is  judged  to  be  essentially  uniform 
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coverage  in  the  impact  zone.  This  position  was  measured  to  a resolution 
of  0.0001  in.  on  the  photomicrograph.  The  impact  locations  were  curve- 
fit  to  an  hyperbola  passing  through  the  sharp  tip  of  the  model  and  the 
model  angle  of  attack  was  established  by  a tangent  to  the  hyperbola  at 
the  nosetip.  In  all  cases,  the  model  angle  of  attack  determined  in  this 
manner  was  greater  than  the  angle  of  attack  measured  before  the  run, 
suggesting  a deflection  due  to  the  air  loads  during  the  run.  Both  angles 
of  attack  are  given  in  Table  3. 

Particle  trajectory  data  for  a typical  case  are  given  in  Fig.  23. 
These  data  are  for  25-P glass  beads  passing  through  a 6.53-deg  Prandtl- 
Meyer  expansion  at  M = 5.  5,  with  pQ  = 996  psi,  HQ  = 437  Btu/lbm.  Based 
on  the  one -dimensional  program  for  this  case,  the  gas  velocity  was  4,302 
ft/ sec  and  the  particle  velocity  was  3,  900  ft/sec.  In  this  case  the  cross- 
stream deflection  of  0.0327  in.  was  slightly  less  than  the  predicted  deflec 
tion  of  0.  0369  in.  at  6.  65  An.  from  the  nose  of  the  model. 

A summary  of  predicted  cross-stream  deflection  6.  65-in.  down- 
stream of  6.  2 5 -deg  Prandtl- Meyer  corner,  for  a large  range  in  particle 


Figure  21.  Photograph  of  Prandtl-Meyer  expansion  test  model  $ 

after  run  in  DET.  | 
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Figure  22.  Photomicrograph  of  impact  bar  from  Prandtl-Meyer  test  model. 
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diameters,  is  given  in  Fig.  24.  The  calculations  predict  all  particles 
with  diameters  greater  than  iO-M  will  deflect  no  more  than  0.  10  in.  , 
representing  nearly  linear  trajectories.  All  particles  having  diameters 
less  than  1 M will  essentially  follow  the  flow  streamlines.  The  transition 
region  covers  a fairly  narrow  decade  of  diameters  between  1 and  10  m. 
The  experimental  points  for  the  beads  and  the  TiC>2  particles  are  in  good 
agreement  with  the  predictions,  but  the  deflection  of  the  MgO  particles 
of  50 -p  diameter  is  about  three  times  the  predicted  value.  Except  for 
the  single  MgO  point,  the  data  establish  the  adequacy  of  the  two-dimen- 
sional trajectory  program  for  design  purposes.  It  is  not  known  whether 
the  MgO  point  is  an  exceptional  point  or  whether  it  reflects  increased 
drag  due  to  the  irregular  shape  of  the  MgO  particles.  In  view  of  the 
good  agreement  of  one -dimensional  acceleration  of  MgO  particles  with 
calculations,  it  is  not  deemed  likely  that  the  discrepancy  is  the  result 
of  irregular  shapes.  The  test  data  appear  to  confirm  the  calculations 
at  the  extremes  of  diameters  near  the  limiting  cases,  but  unfortunately, 
no  particle  samples  were  available  which  would  have  fallen  in  the  tran- 
sition regime. 


Figure  24.  Summary  of  particle  deflection  data  6.65  in.  from 
Prandti-Meyer  corner,  DET,  M = 5.9. 
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4.0  ANALYSIS  AND  DESIGN  OF  A COMBINED 
ABLATION/EROSION  NOZZLE 

An  analysis  of  the  proposed  ablation/ erosion  nozzle  -was  conducted, 
using  as  design  tools  the  programs  for  gas/ particle  flow  described  in 
Section  2.0  and  which  had  been  experimentally  verified  as  described  in 
Section  3.0, 


4.1  PARAMETRIC  DESIGN  STUDY  OF  ABLATION/EROSION 
NOZZLE 


A detailed  numerical  analysis  of  the  nozzle  configuration  shown  in 
Fig.  2 was  made.  The  arc-heater  reservoir  conditions  were  assumed  to 
be  - 200  atm  and  HCJ  = 2,  100  Btu/lbm,  which  are  representative  maxi- 
mum pressure  operating  conditions  for  a typical  high-pressure  arc  heater. 
The  light -gas  particle  acceleration  nozzle  was  assumed  to  intersect  the 
arc-heated  air  nozzle  at  angles  of  10,  15,  and  20  deg.  The  divergence 
half -angle  of  the  light -gas  nozzle  was  assumed  to  be  0.  5 deg.  The  initial 
calculations  were  made  with  helium  as  the  light  gas  and  with  reservoir 
conditions  of  1,000  psi  and  2,0Q0°R. 


4.1.1  Flow-Field  Analysis 

Flow -field  conditions  upstream  and  downstream  of  the  helium  and 
air  shock  waves  were  obtained  by  use  of  the  computer  program  described 
in  Section  2.2,  which  performed  isentrepic  expansion  calculations  up  to 
the  shocks  and  solved  the  oblique  shock  relations  for  the  different  deflec- 
tion angles,  with  the  constraint  that  the  static  pressures  on  both  sides  of 
the  helium/ air  interface  be  equal.  So  that  conditions  downstream  of  the 
air  shock  be  representative  of  a high-pressure  ablation  jet,  it  was  desired 
that  the  pressure  upstream  of  the  air  shock  be  approximately  5 atm 
Parameters  used  as  input  to  the  program  were  the  arc-heater  reservoir 
conditions,  the  air  nozzle  Mach  number,  the  helium  nozzle  reservoir 
conditions,  and  the  impingement  angle  between  the  two  jets.  Results  of 
these  calculations  are  given  in  Fig.  2 5.  Note  that  the  static  pressure 
downstream  of  the  air  shock  is  15.  95  atm  for  the  20 -deg  deflection  case. 
This  impingement  angle  of  20  deg  was  selected  for  further  detailed 
analyses  of  the  configuration,  mainly  because  of  anticipated  design  prob- 
lems with  use  of  smaller  angles. 

4.1.2  Light-Gas  Acceleration  Nozzle  Analysis 


Flow  in  the  helium  nozzle  was  analyzed  with  the  one -dimensional 
particle  acceleration  program  descrioed  in  Section  2.  1,  The  throat 
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Figure  25.  Results  of  flow-field  calculations  for  proposed 
ablation/erosion  nozzle. 

diameter  was  assumed  to  be  0.  119  in.  The  required  static  pressure  in 
the  helium  flow  upstream  of  the  shock  is  0.  82  atm  for  the  20-deg  deflection 
case  (Fig.  25),  and  this  determined  the  length  of  the  helium  nozzle  to  be 
18.2  in.  The  results  of  the  boundary -layer  calculations  indicate  a merging 
of  the  boundary  layer  several  inches  upstream  of  the  18.  2 -in.  length.  In 
view  of  the  uncertainty  of  the  effect  of  the  boundary  layer  merging  on  the 
particle  acceleration,  operation  of  the  nozzle  will  be  attempted  in  this 
configuration.  However,  if  merging  creates  a serious  problem  in  prac- 
tice, it  may  be  necessary  to  increase  the  nozzle  divergence  in  this  region. 
Calculations  were  made  for  MgO  particles  having  various  diameters  be- 
tween 10  and  1,  000  /< , and  the  final  particle  velocities  were  found  to  range 
from  8,  000  to  1,  700  ft/sec.  The  final  gas  velocity  was  10,  150  ft/sec. 

4.1.3  Particle  Trajectory  Calculations 

The  two-dimensional  trajectories  of  particles  crossing  the  four 
flow  field  regions  of  Fig.  25  were  calculated  with  the  program  described 
in  Section  2.  2.  Calculations  were  made  for  the  20-deg  impingement  case 
and  with  MgO  particle  diameters  from  10  to  1,  000  p . As  a conservative 
measure  to  allow  for  possible  boundary -layer  merging,  the  initial  particle 
velocities  for  the  trajectory  calculations  were  assumed  to  be  those  at  the 
n.  ”ge  point  rather  than  at  the  end  of  the  nozzle.  It  was  found  that  there 
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•was  a strong  tendency  for  the  trajectories  downstream  of  the  helium  shock 
to  converge  to  a single  trajectory.  This  results  from  the  fact  that  the  ini- 
tial velocities  developed  in  the  acceleration  nozzle  vary  inversely  with  the 
particle  size  or  mass.  In  Fig.  26,  the  trajectories  of  all  particles  in  the 
range  from  50  to  1,000  p are  shown  as  a single  curve.  The  1 0 -/u  particles, 
however,  showed  substantially  greater  deviation  fron  a straight  line  and 
are  represented  by  a separate  curve.  A summary  of  particle  trajectory 
parameters  is  given  in  Table  4 for  the  five  locations  indicated  in  Fig.  26. 
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Figure  26.  Particle  trajectories  calculated  in  ablation/erosion  nozzle. 


Table  4.  Characteristics  of  Trajectories  in  Ablation/Erosion  Nozzle 


NORMAL 

inches 

OEF. , 

AS, 

degrees 

FARTICLE  VELOCITY 
ft /sec 

LOCATION 

Dp=  10 
b 

100 

b 

1000 

b 

10 

b 

100 

b 

000 

b 

10 

b 

100 

b 

1000 

b 

A- HELIUM  SHOCK 

0 

0 

0 

0 

0 

0 

8028 

4390 

1700 

B-He/AIR  BOUNDARY 

.0014 

.0004 

.0005 

0.62 

0.16 

0.21 

8030 

4405 

1713 

C-AIR  SHOCK 

.0511 

.0158 

.0192 

4.57 

1.34 

1.61 

7918 

4510 

1808 

D-E.I4"  FROM  ~~ 1 

He  SHOCK 

.K34 

.0339 

.0406 

3.74 

1.31 

1.52 

7604 

4547 

1890 

E- 2 86"  FROM 

He  SHOCK 

.1524 

.0505 

.0575 

3.22 

1.28 

1.45 

7427 

4578 
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j It  is  evident  that  particles  having  diameters  from  50  to  1,000  M de- 

flect 0.  05  to  0,  06  in.  normal  to  the  initial  direction  and  the  velocity  vector 
i is  turned  by  only  1.3  to  1.5  deg.  These  trajectories  are  essentially  straight 

lines.  The  nozzle  could  be  used  as  shown  in  Fig.  26  or  the  geometry  could 
j be  tailored  to  the  trajectories  by  a small  rotation  of  the  upper  surface  of 

the  skimmer.  The  linear  and  angular  deflections  of  the  10-M  particles  are 
approximately  three  times  larger  than  for  the  larger  particles.  Since  the 
proposed  facility  concept  will  operate  with  100-m  particles,  this  is  not  a 
handicap  to  the  design. 

4.1.4  Skimmer  Heat  Transfer 

The  skimmer  section  will  no  doubt  present  a difficult  design  problem 
because  of  the  high  heat-transfer  rates.  Calculations  of  the  turbulent 
heating  rates  on  the  air  side  of  the  skimmer  give  values  approaching 
i 10,000  Btu/ft^-sec  - >ar  the  tip,  but  average  values  over  the  first  inch  are 

only  5,700  Btu/ft^-sec.  The  average  heat  load  could  easily  be  absorbed 
by  internal  water  cooling;  however,  the  sharp  tip  of  the  skimmer  may  have 
to  be  rounded  slightly  to  reduce  the  peak  heating  locally.  There  may  also 
| be  some  slight  relief  as  a result  of  laminar  heat  transfer  near  the  tip. 

I 

I 

| 4.1.5  Use  of  Hydrogen  as  Accelerator  Gas 

i The  performance  of  the  acceleration  nozzle  could,  of  course,  be 

increased  by  increasing  the  reservoir  temperature  of  the  helium.  Ultimate 
maximum  performance  of  the  nozzle  could  be  obtained  by  the  use  of  hydro- 
gen, since  it  has  twice  the  specific  heat  of  helium.  Preliminary  calcu- 
lations were  made  for  hydrogen  in  the  present  nozzle,  and  for  the  same 
| light  gas  reservoir  conditions  of  1,000  psi  and  2,000°R,  velocities  of 

| 100-p  MgO  particles  were  indicated  to  be  7,  170  ft/.c;ec,  compared  with 

■ 4,390  ft/sec  in  the  helium  nozzle.  However,  the  boundary -layer  calcu- 

lation indicates  that  a substantially  earlier  merging  and  attainment  of  the 
hydrogen  performance  would  be  contingent  upon  avoiding  any  problems 
associated  with  boundary  layer  merging.  In  addition,  the  use  of  hydrogen 
would  present  many  safety  problems  and  hazards  which  would  require 

i thorough  investigation. 

1 

i 

4.2  DETAIL  DESIGN  OF  A FIRST-GENERATION 
ABLATION/EROSION  NOZZLE 

j 

: For  the  purpose  of  proof -testii .g  a first  generation  ablation/ erosion 

nozzle,  it  was  decided  to  relax  the  arc-heater  operating  conditions  to  150 
atm  in  order  to  avoid  troublesome  operational  problems  with  the  archeater. 
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A design  stagnation  enthalpy  of  2,  200  Btu/  lbm  was  chosen  for  this  pres- 
sure level.  It  was  felt  that  these  conditions  will  provide  a reasonable 
ablation  environment  in  which  to  assess  performance  of  the  nozzle. 
Reservoir  conditions  for  the  light-gas  acceleration  nozzle  were  chosen 
to  be  1,000  psi  and  1,75CPR,  using  helium.  This  also  represents  a 
small  relaxation  in  operating  conditions  from  those  used  in  the  analytical 
feasibility  study  (Section  4.  1).  Velocities  of  100-p  MgO  particles  should 
still  be  in  the  vicinity  of  4,  000  ft/sec,  which  is  considered  adequate  for 
the  proof  test. 

The  design  Mach  number  for  the  arc-heater  nozzle  was  determined 
by  an  isentropic  expansion  from  the  arc-heater  conditions  to  the  desired 
5-atm  pressure  at  the  helium/air  interface.  The  acceleration  nozzle 
Mach  number  was  also  determined  by  the  5-atm  matching  condition  and 
the  helium  reservoir  conditions,  using  the  oblique  shock  relation.  The 
Mach  numbers,  pressures,  temperatures,  and  calculated  particle  veloc- 
ities for  100-p  MgO  are  given  in  Table  : or  the  design  version  ablation/ 

erosion  nozzle. 


Table  5.  Design  Conditions  for  First-Generation  Ablation/Erosion 
Nozzle,  Hot-Flow  and  Cold-Flow  Versions 


GAS 

HELIUM 

Po 

1000  psi 

1750-R 

AIR 

Po 

• 150  atm 

H0 

= P200  B--u. 

lbm 

HELIUM 

Po 

* 73b  psi 

To 

- 500°R 

UPSTREAM  Of  SHOCK 


HOT  po  * '000  psi  3.819  0.818  299  9494  P.058 

FLOW  Tp  - ,750°R 


DOWNSTREAM  OF  SHOCK 


T,  I V,  P ■ 

ft/sec  atm 


4.2.1  Arc-Heater  Nozzle  Design 


A basic  decision  was  made  to  use  two-dimensional  nozzles  in  order 
to  simplify  the  geometry  of  the  intersection  of  the  two  flows.  The  exit  of 
the  arc-heater  nozzle  was  chosen  to  be  0.  70  in.  square.  A nozzle  ex- 
panding only  in  one  direction  was  selected,  giving  a throat  cross  section 
of  0.  189  by  0.  70  in.  for  the  design  Mach  number  of  2.  731.  The  nozzle 
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contour  was  obtained  by  a two-dimensional  method-of-characteristics 
program.  A boundary -layer  correction  was  calculated  but  was  deemed 
to  be  too  small  to  include  in  the  final  contour.  This  nozzle  will  be  cooled 
with  500-psi  demineralized  water  flowing  through  back-side  passages 
(Fig.  27). 


Figure  27.  Cross  section  of  ablation/erosion  nozzle,  hot-flow  version. 


4.2.2  Helium  Particle  Acceleration  Nozzle 

j 

j The  gas  heater,  stilling  chamber,  particle  injection  system,  and 

| particle  acceleration  nozzle  described  in  Section  3.  1.2.  1 were  incorpo- 

I rated  into  the  design  of  the  first-generation  ablation/erosion  nozzle,  with 

j the  exception  that  the  nozzle  was  cut  off  at  the  length  necessary  to  obtain 

j the  correct  static  pressure.  This  length  was  determined  by  calculations 

' with  the  one -dimensional  program  described  previously. 

i 

A round-to-square  transition  section  (0.  364-in.  -diam  to  0.  364-in.  - 
square)  was  designed  to  be  located  between  the  acceleration  nozzle  and  the 
injection  assembly.  Since  this  nozzle  will  be  operated  with  no  cooling, 
provision  must  be  made  for  thermal  expansion.  The  end  of  the  acceleration 
nozzle  is  allowed  to  slide  inside  the  transition  section,  and  while  this  will 
create  some  distrubances  in  the  flow,  it  is  preferable  to  attempting  to 
cool  the  entire  nozzle. 
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4.2.3  Injection  Chamber  or  Assembly 


The  dimensions  of  the  injection  chamber  are  determined  by  the 
exit  dimensions  of  the  arc-heater  nozzle  (0.  7 in.  square)  and  the  helium 
acceleration  nozzle  (0.412  in.  square).  The  mixing  region  between  the 
two  gases  was  originally  assumed  to  be  thin  and  was  neglected.  This 
assumption  is  currently  being  reassessed,  based  on  results  with  the 
cold-flow  nozzle  test  described  in  Section  5.  0.  To  allow  for  uncertainty 
of  the  mixing  layer  thickness  and  to  permit  fabrication  to  proceed,  the 
injection  chamber  was  split  into  two  sections  (Fig.  27).  This  feature 
also  has  the  advantage  of  providing  an  easily  replaceable  piece  in  the 
region  which  is  most  likely  to  receive  particle  erosion  damage. 

Design  of  the  splitter  wedge  to  be  located  in  the  downstream  section 
of  the  chamber  has  not  been  completed.  The  heat-transfer  conditions 
noted  in  Section  4.1.3  suggest  the  use  of  either  a water-cooled  copper 
wedge  or  an  uncooled  wedge  made  of  some  high  temperature  ablation 
material.  In  either  case,  it  will  be  impossible  to  preserve  a perfectly 
sharp  tip  on  the  skimmer  and  an  adjustment  in  position  will  be  required 
to  allow  for  shock  detachment  at  the  nose  of  the  skimmer. 


5.0  DESIGN  AND  TESTING  OF  COLD-FLOW 
ABLATION/EROSION  NOZZLE 


It  was  decided  early  in  the  study  of  the  ablation/ erosion  nozzle 
concept  that  it  would  Je  desirable  to  design,  construct,  and  operate  a 
cold-flow  version  of  this  nozzle  prior  to  final  fabrication  of  the  hot- 
flow  version.  This  would  permit  identification  of  operational  problems 
in  time  to  allow  design  modification  of  the  hot  nozzle,  serve  as  a check 
on  the  design  calculation  methods,  and,  incidentally,  demonstrate  ex- 
perimentally that  particle  transfer  across  the  mixing  layer  was  possible. 

5.1  DESIGN  CONSIDERATIONS 

The  requirement  of  5 -atm  static  pressure  at  the  mixing  layer  be- 
tween the  two  flows  was  retained  in  the  cold -flow  version.  The  throat 
for  the  nozzle  simulating  the  arc-heated  flow  was  chosen  to  have  essen- 
tially the  same  area  as  the  standard  0.  375-in.  -diam  nozzles  used  in 
most  AEDC  arc  heaters,  with  reservoir  conditions  of  735  psi  and  500°R 
to  give  an  air  mass  flow  of  2 lb/sec.  The  design  Mach  number  of  2.  157 
for  this  nozzle  was  determined  by  the  ratio  of  static  to  reservoir  pressure, 
and  this  Mach  determined  the  oblique  shock  angle  on  the  splitter  wedge  for 
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the  20-deg  deflection  angle.  The  particle  acceleration  nozzle  and  particle 
injection  system  were  considered  identical  to  the  hot  nozzle  version.  The 
design  conditions  for  the  cold-flow  version  of  the  ablation/ erosion  nozzle 
are  given  in  Table  5. 


5.2  EXPERIMENTAL  APPARATUS 


ft 

S 


The  components  of  the  cold-flow  nozzle  are  the  cold-flow  air  nozzle, 
the  helium  acceleration  nozzle,  and  the  particle  injection  section  or  cham- 
ber (Fig,  28).  Also  used  in  cold-flow  testing  are  a laser  shadowgraph 
system  and  the  usual  instrumentation  systems.  A photograph  of  the  assem- 
bled system  is  shown  in  Fig.  29. 

5.2.1  Cold- Flow  Air  Nozzle 

The  exit  of  the  cold-flow  air  nozzle  was  0.  5-in.  square,  and  the 
throat  dimensions  were  0.  235  in.  wide  and  0.50  in.  high.  The  contour 
of  this  nozzle  was  obtained  by  a two-dimensional  method -of -characteristics 
program  with  a boundary -layer  correction  added.  This  nozzle  was  fabri- 
cated from  2024 -T4  aluminum,  selected  for  ease  of  fabrication.  It  was 


Figure  28.  Cross  section  of  ablation/erosion  nozzle,  cold  flow  version. 
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Figure  29.  Photograph  of  cold-flow  ablation/erosion  nozzle  apparatus. 

supplied  by  a high-pressure  air  line  or  by  a gaseous -nitrogen  supply  line. 
Either  system  was  capable  of  passing  the  2 -lb/ sec  mass  flow  required  and 
the  stagnation  pressures  were  controllable  from  atmospheric  to  over  750 
psi.  The  design  Mach  number  of  2.  157  is  low  enough  that  liquefaction  was 
not  a concern:  h wevc; , ondensation  of  water  vapor  in  the  air  supply  was 
found  to  be  a problem. 


5.2.2  Acceleration  Nozzle  and  Transition 

It  was  necessary  lorten  the  acceleration  nozzle  to  provide  the 

correct  static  pressur  ■ 94  atm)  at  the  correct  exit  Mach  number  of 

3.43,  As  in  the  hot-flow  version,  a round-to-square  transition  section 
fabricated  from  6061-T6  aluminum  was  provided  between  the  end  of  the 
acceleration  nozzle  and  the  mixing  chamber  or  particle  injection  section. 
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Helium  supply  for  this  nozzle  was  provided  by  a 2,200-psi  helium  trailer, 
Some  testing  was  also  done  with  both  the  air  and  nitrogen  supply  referred 
to  above. 


5.2.3  Particle  Injection  Chamber 


-J 


The  size  of  the  section  in  which  the  particles  were  injected  into  the 
airflow  was  determined  by  the  exit  dimensions  of  the  cold  air  nozzle  and 
the  helium  acceleration  nozzle  (0.  50  and  0.  388  in.  square,  respectively) 
and  by  the  requirement  for  a 20 -deg  impingement  angle  between  the  two 
flows.  Oblique  shock  calculations  yield  a cold  airflow  exit  height  of  0.  321 
in.  The  mismatch  between  this  dimension  and  the  exit  dimension  of  the 
acceleration  nozzle  will  require  that  some  particles  will  impinge  on  the 
upper  wall  of  the  injection  chamber.  Accordingly,  a removable  insert 
was  fabricated  which  allows  for  easy  replacement  of  this  area  (Fig.  28). 

Static  pressure  orifices  were  provided  in  the  injection  chamber 
to  allow  measurement  of  pressure  upstream  and  downstream  of  the 
helium  and  air  shock  waves.  In  addition,  the  sidewalls  of  the  injection 
chamber  were  fabricated  of  removable  sections  of  schlieren-quality 
glass  to  allow  shadowgraph  studies  of  the  flow  field. 

The  stainless  steel  splitter  wedge  was  also  located  m the  injection 
chamber.  Provision  was  made  for  varying  the  angle  of  attack  and  axial 
position  of  the  wedge  to  permit  study  of  the  effect  of  wedge  orientation 
on  the  injection  chamber  flow  field.  The  wedge  position  could  not  be 
varied  during  a run;  however,  rapid  changes  are  possible  between  runs. 


I 

5.2.4  Laser  Shadowgraph  System  j 
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Most  of  the  equipment  used  in  the  attempt  to  obtain  particle  field 
holograms  was  readily  adapted  to  a shadowgraph  system.  For  rapid 
viewing  of  test  results,  a Spectra-Physics  Model  135-5  5-mw  helium- 
neon  laser  was  used  with  high-speed  Polaroid®  black  and  white  film. 

If  these  shadowgrams  were  considered  useful,  an  Apollo  ruby  laser 
was  used  to  obtain  a higher-quality  permanent  record  in  which  consid- 
erably more  flow  details  are  visible. 

5.3  EXPERIMENTAL  RESULTS 

A limited  number  of  shakedown  runs  was  made  with  the  windows 
and  skimmer  wedge  removed  from  the  injection  chamber  assembly. 
Both  helium  and  air  were  used  in  the  particle  acceleration  nozzle  for 
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these  runs.  The  initial  runs  were  made  with  helium  and  resulted  in  a 
slight  erosion  on  the  removable  insert  on  the  top  wall  of  the  particle 
injection  chamber,  apparently  the  result  of  residual  dust  left  in  the 
stilling  chamber  after  previous  runs  of  the  particle  acceleration  nozzle. 
When  the  particle  acceleration  nozzle  was  operated  with  air,  however, 
the  erosion  on  the  insert  was  much  greater,  and  it  was  concluded  that 
this  was  the  result  of  ice  particles  formed  in  the  expansion  of  moist  air. 
The  eroded  area  matched  almost  exactly  the  height  of  the  particle  accel- 
eration nozzle,  thus  giving  inadvertent  proof  of  particle  transfer  across 
the  mixing  layer  in  essentially  linear  trajectories.  In  subsequent  oper- 
ations, the  use  of  moist  air  in  the  light-gas  acceleration  nozzle  was 
avoided. 

For  the  next  series  of  tests,  the  flow  was  properly  confined  to  the 
injection  chamber  by  installation  of  the  windows,  and  runs  were  made 
without  the  skimmer  wedge  to  investigate  the  nature  of  the  mixing  layer 
when  undisturbed  by  the  presence  of  the  wedge.  A shadowgram  of  a 
typical  test  condition  is  given  in  Fig.  30.  The  mixing  region  is  clearly 
seen  and  is  observed  to  propagate  into  the  airflow  region  at  a constant 
rate.  Little  penetration  into  the  helium  flow  from  the  acceleration  nozzle 
is  noted.  Some  evidence  of  separation  is  observed  in  the  airflow  upstream 
of  the  junction  with  the  helium  flow  and  a weak  oblique  shock  results.  The 
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mixing  region  has  a considerable  thickness  even  at  the  far  upstream  end 
of  the  region.  It  was  also  observed  during  these  tests  that  the  pressure 
measured  at  the  upstream  pressure  tap  in  the  acceleration  nozzle  did 
not  match  the  calculated  value  and,  in  fact,  closely  matched  the  value 
calculated  for  the  region  downstream  of  the  shock.  This  indicated  that 
the  oblique  shock  in  the  helium  flow  was  displaced  upstream  of  its  cal- 
culated position. 

The  skimmer  wedge  was  next  installed  at  its  intended  location,  and 
the  cold-flow  ablation/erosion  nozzle  was  operated  with  nitrogen  and  helium 
at  near  design  conditions.  A shadowgram  of  the  flow  with  the  wedge  in  this 
position  is  given  in  Fig.  31.  A prominent  mixing  layer  is  noted,  with  a 


Figure  31.  Shadowgraphs  of  flow  in  cold-flow  nozzle,  with  wedge 

skimmer,  normal  wedge  position,  and  adjusted  wedge  position. 
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relatively  strong  oblique  shock  originating  at  the  beginning  of  the  mixing 
layer.  The  mixing  layer  is  deflected  into  the  main  airflow  stream,  which 
apparently  narrows  the  flow  channel  and  causes  subsonic  flow  in  this  region. 
It  was  assumed  that  proper  adjustment  of  the  operating  pressures  of  the 
two  nozzles  would  restore  the  intended  flow  configuration. 


The  airflow  nozzle  stagnation  pressure  was  held  near  the  design 
condition,  and  the  helium  nozzle  stagnation  pressure  was  lowered  to  400 
psi.  Although  the  upstream  static  pressure  in  the  airflow  did  not  indicate 
that  supersonic  flow  had  been  established,  examination  of  shadowgrams 
revealed  a thinner  mixing  layer  and  a more  regular  shock  structure.  The 
mixing  layer,  however,  was  still  deflected  ahove  the  wedge.  Analysis  of 
both  the  pressure  measurements  and  the  shadowgrams  suggested  that  the 
wedge  was  acting  like  a two-dimensional  blockage  in  a supersonic  wind 
tunnel.  It  was  decided  to  move  the  wedge  downstream  and  tilt  it  upward 
until  the  area  required  to  allow  supersonic  flow  above  the  wedge  was 
achieved.  In  addition,  some  material  was  removed  from  the  lower  surface 
of  the  wedge  in  an  effort  to  improve  the  flow  in  the  helium  leg  downstream 
of  the  wedge  tip.  A shadowgram  of  the  flow  made  after  these  changes  shows 
considerable  improvement  in  the  flow  field  (Fig.  31).  The  mixing  region 
is  once  again  the  thickness  attained  during  the  tests  without  the  wedge,  and 
it  now  intersects  the  wedge  at  the  proper  location.  The  shock  structure  in 
the  airflow  region  of  the  injection  chamber  is  close  to  design,  although  it 
is  located  too  far  downstream  to  intersect  the  corner  of  the  insert  on  the 
upper  surface.  The  Mach  number  calculated  from  the  20 -deg  deflection 
angle  and  the  shock  angle  on  the  wedge  is  2.  048,  compared  with  the  design 
value  of  2.  157.  This  reduction  in  Mach  number  was  presumed  to  be  the 
result  of  blockage  by  the  mixing  layer.  An  effort  was  made  to  calculate 
this  blockage  effect  with  an  existing  mixing  layer  program,  and  a Mach 
number  of  2.05  was  the  result,  very  close  to  the  measured  value.  (Note 
that  this  reduction  in  Mach  number  is  not  the  result  of  boundary -layer 
displacement  in  the  air  nozzle,  since  a boundary -layer  correction  was 
added  to  the  inviscid  nozzle  contour. ) 

The  static  pressures  measured  at  the  "upstream"  orifice  in  the  I 

helium  flow  continued  to  indicate  that  the  shock  wave  was  upstream  of  ! 

this  location.  If  it  is  assumed  that  the  shock  position  is  just  upstream  j 

of  this  orifice,  it  is  possible  to  calculate  new  values  of  Mach  number 
and  static  pressure  in  the  helium  flow  from  the  oblique  shock  equations.  1 

These  values  match,  within  a few  percent,  values  obtained  during  oper-  j 

ation  of  the  light -gas  acceleration  nozzle  alone.  This  indicates  that  the  j 

transition  section  is  functioning  properly  and  that  the  oblique  shock  is  at  j 
an  angle  greater  than  design  for  some  reason. 
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To  benefit  from  the  experience  trained  in  operation  of  the  cold-wall 
version  of  the  ablation/erosion  nozzle,  the  design  of  the  hot-flow  version 
will  be  slightly  modified  to  allow  for  the  greater  than  expected  blockage 
of  the  mixing  layer.  The  hot-flow  version,  from  the  beginning,  has  had 
the  requirement  that  the  dimensions  of  the  final  ablation/erosion  jet  be 
equal  to  the  exit  of  the  light-gas  acceleration  nozzle,  to  minimize 
particle  erosion  within  the  nozzle. 

6.0  CONCLUSIONS 

Results  of  a test  to  analyze,  design,  and  test  the  components  of  a 
combined  ablation/ erosion  nozzle  may  be  summarized  as  follows: 

1.  An  analytical  study  was  made  to  determine  the  possibility 
of  attaining  a combined  ablation  and  erosion  test  environ- 
ment in  a facility  which  has  operating  times  measured  in 
seconds. 

2.  Computer  programs  were  developed  and/or  extended  for 
analysis  of  one-dirnensional  particle/gas  flows  with  thick 
boundary  layers  and  for  analysis  of  two-dimensional  particle 
trajectories  in  various  gas  flow  fields.  An  empirical  particle 
drag  coefficient  relation  was  developed  based  on  most  recent 
test  data. 

3.  It  was  determined  that  direct  drag-acceleration  of  100-P  -diam 
particles  in  typical  low  Mach  number  nozzles  used  in  arc- 
heated  ablation  facilities  was  relatively  ineffective.  Such 
particles  can  be  accelerated  to  hypervelocities  in  a long, 
slowly  diverging  nozzle  in  which  length  of  exposure  to  high 
drag  force  is  maximized. 

4.  By  adopting  the  concept  of  acceleration  of  particles  in  a 
separate  light-gas  acceleration  nozzle,  followed  by  injec- 
tion into  an  ablation  jet  by  inertial  transfer  across  a mixing 
layer  between  the  two  jets,  detailed  analysis  indicates  an 
erosion  capability  can  be  added  to  a state-of-the-art  ablation 
jet. 

5.  A detail  design  was  completed  for  a facility  based  on  this 
concept,  and  a cold -flow  version  was  built  and  operated  to 
obtain  operating  experiem  '■  with  the  new  configuration. 

Particle  velocities  which  were  measured  in  the  acceleration 
nozzle  agreed  well  with  the  predictions  of  the  one -dimensional 
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program.  Aerodynamic  calibration  of  the  combined  flow  indi- 
cated that  the  displacement  thickness  of  the  mixing  layer  was 
substantial  and  should  be  allowed  for  in  final  design. 
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A 

ap 

acv  al»  a2> 
bo,  bi»  b2» 
co»  ci,  c2, 
d0,  d1#  d2 

CD 

Cf 


Cl,  C2,  C3 
D 

Ho 

h 

M 

mp 

Nu 

Po 

Pr 

P 

Qp 

Re 


NOMENCLATURE 


Flow  area 

Cross-sectional  area  of  particle 

Coefficients  in  drag  coefficient 
correlation,  Eq.  (2)  and  Table  2 

Particle  drag  coefficient 
Skin  friction  coefficient 
Stanton  number 
Specific  heat 

Constants  in  Eq.  (5)  for  skin  friction 
Diameter 

Total  enthalpy  based  on  energy  balance 
Static  enthalpy 

Mach  number  of  flow  or  relative  flow 
Mass  of  particle 
Nusselt  number 


I 

I 

* ! 


1 

» 1 


1 

i 

1 

\ 

I 

i 


I 


3 

.•1 

\ 

4 


Impact  pressure 
Prandtl  number 
Static  pressure  of  gas 
Heat  flow  to  particle 
Reynolds  number 
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V 

■A 

■1 
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i i! 

;! 


K « 


SP 

Particle  surface  area 

T 

Temperature 

V 

Velocity 

V 

Volume  of  particle 

X 

Streamwise  coordinate 

a 

Angle  of  attack  or  angle  of  Prandtl-Meyer 
expansion  ray 

y 

Specific  heat  ratio 

6 

Total  thickness  of  boundary  layer 

6* 

Displacement  thickness  of  boundary  layer 

1 

$ 

j 

Momentum  thickness  of  boundary  layer  or 

angle  of  velocity  vector  \ 

P 

Viscosity 

P 

Dens  ity 

TT 

Barred  quantity  is  an  average  value  across 
the  flow 

( ) 

Unbarred  quantity  is  value  in  the  isentropic 
core  of  the  flow 

r> 

Vector  quantity  j 

i 

| 

i 

SUBSCRIPTS  I 

calc 

Calculated  value  j 

inc 

Incompressible  value  ! 

meas 

Measured  value  | 

*4 

o 

Reservoir  condition  I 
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P 

t 

W 


Particle 
Total  condition 
Wall  condition 


